Effects of disorder on the wave front depinning transition in spatially discrete systems 
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Pinning and depinning of wave fronts are ubiquitous features of spatially discrete systems de- 
scribing a host of phenomena in physics, biology, etc. A large class of discrete systems is described 
by overdamped chains of nonlinear oscillators with nearest-neighbor coupling and subject to random 
external forces. The presence of weak randomness shrinks the pinning interval and it changes the 
critical exponent of the wave front depinning transition from 1/2 to 3/2. This effect is derived by 
means of a recent asymptotic theory of the depinning transition, extended to discrete drift-diffusion 
models of transport in semiconductor superlattices and confirmed by numerical calculations. 
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Phenomena in many different fields may be described 
by means of spatially discrete systems: motion of dislo- 
cations in crystals ^ , atoms adsorbed on a periodic sub- 
strate 0, arrays of coupled diode resonators ||], weakly 
coupled semiconductor superlattices (SL) , sliding of 
charge density waves (CDW) superconductor Joseph- 
son array junctions Q , propagation of nerve impulses 
along myelinated fibers , pulse propagation through 
cardiac cells , calcium release waves in living cells Q , 
etc. In many of these systems, disorder due to differences 
in the parameters of individual elements is important be- 
cause it has a strong impact in the collective behavior. A 
distinctive example of collective behavior in discrete sys- 
tems (not shared by continuous ones) is the phenomenon 
of wave front pinning: for values of a control parame- 
ter in a certain interval, wave fronts joining two different 
constant states fail to propagate When the control 
parameter surpasses a threshold, the wave front dcpins 
and starts moving The existence of such thresholds 
is an intrinsically discrete fact, which is lost in contin- 
uum aproximations. Recently, a theory of front depin- 
ning and motion near threshold has been proposed by 
some of us for one-dimensional (ID) nonlinear spatially 
discrete reaction-diffusion (RD) systems In our the- 
ory, propagation failure and front depinning are charac- 
terized by studying the behavior of a few sites, provided 
the effects of spatial discretization are sufficiently strong 

10. 

In this paper, we consider the effect of weak disorder 
on the wave front depinning transition in spatially dis- 
crete ID systems. Applications will include discrete RD 
systems subject to a random field, sliding CDW and do- 
main motion in SL. In these examples, the main effect of 
disorder is to soften the transition, changing the critical 
exponent from 1/2 to 3/2. The latter value was obtained 
by D. Fisher in a mean field model of sliding CDW using 
scaling arguments ]T^ . 



We consider chains of diffusively coupled overdamped 
oscillators in a potential V, subject to a random force 
fieldF-|-7^„: 

du 

— ^ = Un+l - 2u„ + Un-1 + F - Ag{Un) + 7Cn- (1) 

at 

Here g{u) — V'{u) presents a 'cubic' nonlinearity, such 
that Ag{u) - F has three zeros, Ui{F/A) < U2{F/A) < 
U^{F/A) in a certain force interval {g' {Ui{F / A)) > for 
i = 1,3, g'{U2{F/A)) < 0). The fluctuating part of the 
force field is 7<^„, where 7 > characterizes the disorder 
strength and ^„ is a zero mean random variable taking 
values on an interval (—1, 1) with equal probability. An 
example of a model described by Eq. (nl) is (except for 
the mean field approximation which we do not make) 
D. Fisher's modification of the Fukuyama-Lee model of 
sHding CDW 0. In it, u„ = 6'„ — Xn, g{u) = sinu, 
761 = Xn+i - 2xn + Xn-1, whcrc 0n is the CDW phase 
at the site n and Xn is a random variable taking values 
with equal probability on (0, 27r). 

Provided g{u) is odd with respect to C/2(0) and 7 = 0, 
there is a symmetric interval \F\ < F^ where the wave 
fronts joining the stable zeros Ui{F/A) and Uy,{F/A) 
are pinned. For > Fc, there are smooth traveling 
wave fronts, Un(t) = u{n — ct), with u{—oo) = Ui and 
m(oo) — U3. The velocity c{A, F) depends on A and F 
and it satisfies cF < and |c| cx (|F| — Fc) 2 as |F| — > Fc 
|lTt . Examples are the overdamped Frenkel-Kontorova 
(FK) model {g = sinu) ||l^ and the quartic double well 
potential {V = (m^ — 1)^/4). Less symmetric nonlin- 
earities yield a non-symmetric pinning interval and our 
analysis of the depinning transition applies to them with 
trivial modifications 

Let us recall the main features of the active point the- 
ory of the wave front depinning transition in the absence 
of disorder [^,^. Except when A is too small (the con- 
tinuum limit in which Fc 0), the stable wave front 
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profile differs appreciably from either Ui or U3 at finitely 
many points, n = —L, . . . , — 1, 0, 1, . . . , M, called the 
active points. At n < — L, m„ w Ui{F/A) and at n > Af, 
Un ~ J73(F/A). We shall reconstruct the wave front pro- 
file u{n~ct) by analyzing the behavior of the active points 
Un{t), n = —L, . . . , — 1, 0, 1, . . . , M, as the front moves for 
F > Fc > (the case F < can be obtained by using 
symmetry considerations). The arbitrary phase of the 
(translation invariant) wave front will be fixed by impos- 
ing that the solution of the system of active points at 
time t = (and F slightly larger than F^) be equal to its 
stationary solution at F = Fc, Un{A, Fc), n = —L, . . . ,M 
up to terms of order (F — Fc). Fc is obtained from the 
condition that the matrix of the coefficients in the system 
of active points linearized about the stationary solution 
has a zero eigenvalue. Provided Vn (K?^ -!-...-(- V^j — 
1) is the corresponding eigenvector, an outer approxi- 
mation to the solution Un{t) is Un{t) Un{A, Fc) + 
v{i)^n, where the amplitude ip obeys the equation 
d(p/dt = a{F - Fc) + (3ip^, in which a = Ef=-L ^^ + 
A-^[V-L/g'{UiiFc/A)) + VM/g'iUsiFc/A))] > 0, /? = 
-i^/^)J2fL-L9"iMAFc))V^ > 0. The solution 
of this equation such that i^(0) = [equivalent to 
M„(0) = UnjA Fc) for n ^ ~L,..., M] is tp ^ [a{F ~ 
Fc)/P] 5 tan (^a/3(F - Fc)t) . The amplitude (p blows up 
at times t — ±tf,, tb — Tr/{2^/aP{F — Fc)). The inverse 
width of this time interval yields an approximation for 
the wave front velocity, |c| '-^ y/a(3{F — Fc)/tt. At the 
blow up times, the previous outer approximation to u„(t) 
has to be matched to an appropriate inner solution. At 
the later blow up time the appropriate inner solution 
is the solution of the active point system at F = Fc with 
the boundary conditions that u„ — it„(A, Fc) est —00 
and Un = Un+i{A, Fc) as t — > 00. At the earlier blow up 
time —tb, the inner solution obeys the same system of 
equations at F = Fc, but the boundary conditions are 
Un = Un-i{A,Fc) as i — > —CO and w„ — Un{A,Fc) as 

Effects of disorder. How does weak disorder modify 
this picture of the wave front depinning transition? Our 
main idea is to find a dominant balance of the disor- 
der effects with nonlinearities and (F — Fc) near the de- 
pinning transition. Given our active point construction 
sketched above, the dominant balance is struck provided 
(F — Fc) = 0(7) as 7 ^ 0. The amphtude Equation 
becomes 
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and the matching condition is the same as before. The 
solution of Eq. (g) blows up at the end of time intervals 
of duration l/|c/f |, where 



provided the argument of the square root is positive. 
Otherwise the motion of the wave front stops and it be- 
comes pinned. Notice that we have chosen now Ufl(t) as 
the central active point that was distinguished with the 
subscript zero in our previous formulas. After the blow 
up, Ufl jumps to uji^i{A,Fc), approximately, and it re- 
mains there until a time l/jcfl+il has elapsed. Then it 
jumps to uii.^2{A^Fc) approximately, and so on. 

The magnitude of interest in these systems is usually 
an average velocity \cji\ over sufficiently many points. 
For example, this magnitude is proportional to the cur- 
rent due to a sliding CDW and it is important to know 
its behavior near the depinning field and the magnitude 
thereof. We shall argue that the average velocity is ap- 
proximately given by the following equation: 



N 



^Ei^«i = (i^«(^)i)' (4) 



R=l 



(|ci?,(OI> = TT / - + (5) 



Here N ^ {L + M + 1) is sufficiently large, a — 1, and 

{x} '^ is ^/x if a; > and zero otherwise. The idea of 
the proof is as follows. Let us assume that A is so large 
that L = M = and there is only one active point. 
Then the central limit theorem applied to Eq. (^) with 
Vr = ^ yields Eqs. (^ - (||). Let us assume now that 
there are two active points. The previous argument fails 
because now the sum in Eq. (|^) comprises two terms in- 
stead of one. Then the terms in the arithmetic mean are 
no longer independent: when uji — ui for instance, Eq. 
(^ contains ^1 and ^2- After the blow up time, un = U2 
and Eq. ^ contains ^2 and ^3, etc. However, we can 
group the sums appearing in the arithmetic average of 
Eq. (|4|) in two groups containing only independent ran- 
dom variables: R — 2r—l and R = 2r, with r = 1, 2, . . .. 
The variable Vi^i -I- V2^2 has zero mean and correlation 
2cr^/3, where = + V2 = L This correlation is 
exactly the same as that of the variable ^1. Then the 
central limit theorem applied to each group (of sums of 
'dimer' random variables) gives one half the integral in 
Eq. (H), and the sum of these two halves yields Eq. (||). If 
we have more active points, we just have to subdivide the 
arithmetic mean in as many subgroups as active points 
and use the previous argument to prove Eq. (^). 

The elementary integral in Eq. (^ yields c if 
F < Fc — 7cr/a (recall that a — 1), 
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■^i\F-Fc\ < ^, 
~|f(F-Fc)-7|* 



(6) 



\cr\ 



M 



^/3(F - Fc) + 7/5 E VR+nU+n, (3) 



if (F — Fc) > ja/a. Clearly we have a new critical field 
F* — Fc~^a/a, and a new critical exponent 3/2 (instead 
of 1/2 for the case without disorder), |c| cx (F — F*)-^^^. 
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We have compared our theory to the direct numerical so- 
lution of Eq. d) in Fig. H obtaining an excellent agree- 
ment of theoretical predictions and numerical simulation. 

Similarly, we can analyze the effect of weak disorder in 
the doping of the wells on the motion of wave fronts in dc 
current biased semiconductor SL. If the total current den- 
sity is close to a pinning value, the displacement current 
is almost zero except at certain times during wave front 
motion at which the wave front jumps from one well to 
the adjacent one. The average velocity given by Eq. is 
proportional to the arithmetic mean of time averages of 
the displacement current over the time interval between 
maxima thereof. The average velocity is basically the 
mean velocity at which the wave front traverses wells. 
To calculate it, we take advantage of our theory of wave 
front motion in current biased SL The equations 
we use are those in Ref. Q except that the dimensionless 
Poisson equation is now Ei — Ei-i = v(ni — l — "f(,i) where 
7 and ^ defined as in Eq. ([^) represent the disorder in well 
doping. In the SL equations, the roles of the force F and 
the parameter A are taken by the total current density J 
and the dimensionless doping u. If 7 = and v surpasses 
a certain minimal value, there are two critical values of 
the current, Ji and J2, such that a wave front is pinned 
if Ji < J < J2 and it moves at a constant velocity c( J, 
otherwise. The velocity c is positive if J < Ji and neg- 
ative if J > J2. Furthermore, near the critical currents, 
|c| oc I J — Jc|^/^ (Jc is either Ji or J2) (s). How does the 
disorder correct this picture? The effect of disorder is to 
add a term jD{Ej)^j+i--f[v{Ej) + D{Ej)]£,j to the total 
current density J in the dimensionless discrete Ampere 
equation. Then the theory in Ref. yields an equation 
similar to Eq. (ph for the front velocity: 



\cr\ = - . 
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^(3{J - J,) - 7/3 E U],+nWR+n. (7) 
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Wli+n — (w_R,+„ + -DiJ+n) C_R+n — -D/J+nC-R-t-n-t-l • (8) 

Here: (i) C/jj^„ is the left eigenvector corresponding to 
the zero eigenvalue of the linearized equations about 
the stationary solution at J = (chosen so that 
Y.n=-L UR+nUR+n = 1; U R+„ is the right eigenvector), 
(ii) we define VR+n — v{EB.+n), etc., and (iii) a and (3 
are given by Eqs. (20) and (21) in Ref. |||. The val- 
ues Eji+n are those of the stationary electric field profile 
at J = Jc. The noise term in Eq. (^ can be written 

as Y.[UR+n{'"R+n + DR+n) - U]^^^_^D R+n-A^R+n, pro- 
vided we take C/fl_|.„ = for ti < — L and n > AI + 1. 
The noise term has zero average and a correlation 2cr^/3, 

with a2 = EiUj^+JvR+n + DR+n) " f/^+„_ii?i?+„_l]2. 

Using the previously mentioned argument of splitting the 
arithmetic mean in groups of independent random vari- 
ables, we can show that Eqs. (|^) and (^) hold provided a 
in Eq. (|^) is given by the previous formula and {F — Fc) 
is replaced by (J — Jc). Comparisons of the resulting 
average front velocity c with the results of numerically 



solving the SL model are shown in Figures I and I for 
currents close to Ji and J2, respectively. 

In conclusion, we have shown that weak disorder 
changes qualitatively the wave front depinning transition 
in overdamped one-dimensional discrete models. Disor- 
der shrinks the pinning interval and it changes the critical 
exponent for the velocity from 1/2 to 3/2. Whether these 
features are robust and hold for strong disorder remains 
to be seen. An interesting indication is that the criti- 
cal exponent 3/2 is obtained independently of the noise 
strength in mean field models of sliding CDW . 
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FIG. 1. (a) Average velocity |c| as a function of F for 
A = 10, Fc = 6.102281 and 7 = 0.1. (b) Graph of ln|c|/lnF 
showing the crossover to the critical exponent 3/2. 
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FIG. 2. Dimensionless average wave front velocity for the 
9/4 SL with dimensionless parameters f = 3, Ji = 0.179203, 
7 = 0.01. 
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FIG. 3. Same as Fig, 
J2 = 0.790203, 7 = 0.01. 
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